"0","species <- c(LETTERS[1:25])
probs <- c(.9,rep(((1-.9)/24),24))
s.probs <- data.frame(species,probs)

sampled.inds <- table(sample(s.probs$species,100,prob = s.probs$probs, replace =TRUE))
Rarefied.richness <- ddply(inds,.(N),each(EsN.calc))
lowE.highS <- Rarefied.richness
lowE.highS$model <- rep(\"low,high S\",nrow(lowE.highS))

all.sims <- rbind(lowE,intermediateE,highE,highE.lowS,lowE.highS)

limits <- aes(ymax = EsN + Es.se,ymin=EsN - Es.se)
p <- ggplot(all.sims,aes(N,EsN,colour=model))
p  + geom_line() + geom_point() + geom_errorbar(limits,width=0)"
